arXiv:astro-ph/9701094vl 15 Jan 1997 


A&A manuscript no. 

(will be inserted by hand later) 


Your thesaurus codes are: 

06(08.02.1, 08.02.2, 08.14.2, 08.13.1, 08.09.2 HU Aqr (=RXJ2107.9 0518), 

02 . 01 . 2 ) 


ASTRONOMY 

AND 

ASTROPHYSICS 

1.2.2008 


Phase-resolved high-resolution spectrophotometry of the 
eclipsing polar HU Aquarii 

Axel D. Schwope 1 ’ 0, Karl Heinz Mantel 2 ’* *, and Keith Horne 3 

1 Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 14482 Potsdam, FRG (e-mail: ASchwope@aip.de) 

2 Universitatssternwarte Miinchen, Scheinerstrasse 1, 81679 Miinchen, FRG 

3 Univ. of St. Andrews, School of Physics and Astronomy, North Haugh, St. Andrews, Fife KY16 9SS, Scotland, UK 
Received, accepted 


Abstract. We present phase-resolved spectroscopy of the 
bright, eclipsing polar HU Aqr obtained with high time 
(~30sec) and spectral (1.6A) resolution when the system 
was in a high accretion state. The trailed spectrograms 
reveal clearly the presence of three different line compo¬ 
nents with different width and radial velocity variation. 
By means of Doppler tomography their origin could be 
located unequivocally (a) on the secondary star, (b) the 
ballistic part of the accretion stream (horizontal stream), 
and (c) the magnetically funnelled part of the stream (ver¬ 
tical stream). For the first time we were able to derive a 
(near-)complete map of the stream in a polar. We pro¬ 
pose to use Doppler tomography of AM Herculis stars 
as a new tool for the mass determination of these bina¬ 
ries. This method, however, still needs to be calibrated by 
an independent method. The asymmetric light curve of 
the narrow emission line originating on the mass-donating 
companion star reveals evidence for significant shielding of 
60% of the leading hemisphere by the gas between the two 
stars. 


Key words: Accretion cataclysmic variables - AM 
Herculis binaries - stars: HU Aqr (= RX J2107.9-0518) - 
stars: eclipsing 

1. Introduction 

Polars or AM Herculis binaries are strongly magnetic cata¬ 
clysmic binaries. The white dwarf’s magnetic field locks 
both stars in synchronous rotation and prevents the for¬ 
mation of an accretion disk. Instead mass transfer occurs 
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via an accretion stream which is thought to follow initially 
a ballistic trajectory and later governed by the magnetic 
field, which leads the matter along magnetic field lines to 
the polar regions of the white dwarf. Most of the gravi¬ 
tational energy is released there, although some dissipa¬ 
tive energy loss is expected to occur all along the stream, 
in particular in the coupling region where the stream is 
forced to leave the orbital plane. 

The stream manifests itself observationally in mainly 
three ways, (A) by so-called absorption dips observed in 
IR, optical or X-ray light curves (provided the inclination 
is high enough and the observer and the accretion spot 
are (blue in color) located on the same side of the or¬ 
bital plane), (B) by residual emission during eclipses of 
the white dwarf, and (C) by strong, asymmetric and vari¬ 
able emission lines (H, Hei, Hen and other high-ionized 
species) at optical and UV wavelengths (provided the sys¬ 
tem is in a high accretion state). The dips mentioned in 

(A) occur, if the line-of-sight crosses the out-of-plane parts 
of the stream. They are caused by photoelectric (at X-ray 
wavelengths), grey or free-free (at optical/IR wavelengths) 
absorption and are richly structured, indicating a highly 
fragmented stream. An overview of relevant data and anal¬ 
ysis was given recently by Watson (1995). Phenomenon 

(B) is less studied because only a handful of eclipsing po¬ 
lars are known and the numbers of photons collected at 
eclipse which certainly arise from the stream are rather 
few, hence, modelling of these events is widely missing 
in the literature, but nice observational data meanwhile 
have been obtained e.g. for UZ For (Schmidt et al. 1993), 
WW Hor (Beuermann et al. 1990) and HU Aqr (Schwope 
et al. 1995). The emission lines (C) have attracted by far 
most attention because their occurence is one of the ‘clas¬ 
sical’ identification criterion. A diverse appearance is re¬ 
ported in the literature, they may consist of broad, narrow, 
high- and medium-velocity as well as (quasi-)stationary 
components, all with different radial velocity amplitudes 
and systemic velocities. These components may have a 
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different appearance when a specific system is re-observed 
occasionally (see e.g. the review of Mukai 1988). The mini¬ 
mum distinction usually made discerns between a narrow 
line attributed to the secondary star and a broad underly¬ 
ing component whose origin is assumed to be somewhere 
in the accretion stream (see e.g. Liebert & Stockman 1985, 
who have shown that for the systems known at that time 
the radial velocity amplitude of the narrow line lies be¬ 
tween that expected for the L\ and the center of mass 
of the companion). Later on, Rosen et al. (1987) have 
presented high-resolution spectra of V834Cen and distin¬ 
guished 4 subcomponents, which were thought to origi¬ 
nate at different parts of the accretion stream and the 
secondary star. Schwope (1991) and Schwope et al. (1991, 
1993a) made use of repeated measurements of the narrow 
line in order to derive long-term ephemericles of the sec¬ 
ondary star for the non-eclipsing systems QQ Vul, MRSer 
and V834 Cen. Recent attempts to locate the line emission 
in polars applied tomographic techniques to data obtained 
on VV Pup and RXJ0515.6+0105 (Diaz & Steiner 1994, 
Shatter et al. 1995). 

HU Aqr is a system which has all of the observable 
features of the stream mentioned, at least occasionally. It 
was discovered independently by British and German as¬ 
tronomers during their optical identification programmes 
of ROSAT WFC and PSPC sources detected in the cor¬ 
responding sky surveys (Schwope et al. 1993b, Hakala et 
al. 1993). Since it was the brightest eclipsing polar with 
the most extended eclipse at that time (meanwhile sur¬ 
passed by RX J0515.6+0105) it attracted immediate inter¬ 
est. We initiated an intensive optical follow-up study using 
4 telescopes simultaneously and performed spectroscopy 
with high and low spectral resolution, high-speed photom¬ 
etry simultaneously in UBVRI-colours as well as standard 
one-channel photometry using a V-filter. All data have 
good signal/noise-characteristics and provide a compre¬ 
hensive database suitable to reach deeper insight into the 
kinematics and accretion processes of the whole class. Here 
we present as a first paper in an upcoming row our results 
of high-resolution spectroscopy. 

2. Observations and reductions 

HU Aqr was observed with the 3.5m-telescope and 
Cassegrain double-beam spectrograph (TWIN) at Calar 
Alto, Spain, on August 17/18, 1993. Both arms of this 
double beam spectrograph (beamsplitter at ~5500A) were 
equipped with low-noise Tek-CCDs (IK, pixel size 24 urn) 
as detectors. We used a l'/8 wide slit thus matching the 
approximate size of the seeing disk. The spectral cover¬ 
age was 418CU5070A in the blue arm, 7450-8770A in the 
red arm. The red spectra were taken in order to search 
for the photospheric Na-doublet from the secondary star. 
Since HU Aqr was in a high accretion state at that time 
this attempt failed. We, therefore, restrict our analysis to 
the spectra taken in the blue arm. The grating used there 


with 1200 grooves per mm yielded a FWHM spectral res¬ 
olution of ~1.6A in the center, and ~3A at the edges of 
the spectra. 
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Fig. 1. (upper panel) Phase-averaged optical B-band 
lightcurve of HU Aqr for the period August 16-18, 1993. The 
counts of HU Aqr were reduced to that of a simultaneously 
observed nearby comparison star. Binsize of the lightcurve is 
3.75 sec (2000 phase bins), (lower panel) Phase coverage of the 
individual trailed spectra of length 30 min (no. 1-6) and 12.5 
min (no. 7). 


In order to avoid velocity smearing a radial velocity 
study of an AM Herculis binary requires a phase resolu¬ 
tion of A0 ~ 0.05 (Mukai 1988), hence ^6 m for a binary 
like HU Aqr with period P or b = 125 m . The time over¬ 
head between subsequent exposures at the TWIN, on the 
other hand, is ~5 m due to CCD-readout, data transfer 
and storage of the large-format images from both arms of 
the spectrograph, which prevents a dense phase coverage. 
We resolved the time conflict by trailing the object along 
the slit while exposing. Individual exposure times were 
30 m with a nominal trailing rate of 15 sec/pixel, which, 
taking into account the size of the seeing disk of ~2 pixels, 
finally yielded a time resolution of ~30 sec or 0.004 phase 
units. A total of 6 spectra with 1800 s and one with 750 s 
integration time were obtained yielding full-phase cover¬ 
age. Wavelength calibration of the flatfielded and bias- 
subtracted two-dimensional images was performed using 
Helium-Argon arc spectra taken every hour, a small shift 
of the arc line pattern on the CCD due to flexure of the 
spectrograph was accounted for during reduction of the 
data. Cosnric-ray removal was achieved by median filter- 
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ing small regions of size 5x5 pixels around suspected 
cosmics which were identified by Laplace-filtering of the 
raw images. Sky-subtraction was performed afterwards in 
the usual manner by fitting polynomials to CCD-columns 
above and below the region exposed to the target star. 

Next, the spatial coordinate of individual pixels along 
the slit had to be transformed to time or orbital phase. We 
defined start and end of the exposure as coniciding with 
those pixels lying most nearly to half intensity between av¬ 
erage sky- and (sky+object)-level in an adjacent 10 pixel 
surrounding. This resulted in a trailed spectrum with 116 
or 117 pixels for an 1800 sec exposure, slightly different 
from the 120 pixels expected from the nominal trailing 
rate of the telescope. Since the lightcurves extracted from 
the trailed spectrograms and the simultaneous photome¬ 
try yielded no systematic offset, we assigned the discrep¬ 
ancy to a slightly incorrect trailing rate of the telescope 
and estimate the timing accuracy of the spectrum in a 
single CCD-row to be better than 15 sec. 

Simultaneously to our spectroscopy we obtained high¬ 
speed photometric data (0.5 sec resolution) using the 
2.2m-telescope also located at Calar Alto which was 
equipped with the MCCP (Barwig et al. 1987). The 
MCCP allows monitoring of the target star, a comparison 
star and the nearby sky simultaneously in UBVRI colours. 
In order to correct slit-losses and guiding errors of our 
trailed spectra, we first flux-calibrated them using obser¬ 
vations of the white-dwarf standard star Wolf 1346, folded 
them through the response of the MCCP B-filter, binned 
the MCCP-data in 15 sec time bins according to the tim¬ 
ing of our trailed spectra and finally multiplied individual 
rows of the trailed spectra using the ratio of both light 
curves (photometric, spectroscopic). 

Final steps of data reduction included normalization 
of the spectra by subtracting low-order polynomial fits to 
the continuum and phase-folding into 100 phase bins using 
our updated photometric ephemeris of eclipse center 

Ted = HJD2448896.543707(27) + E x 0.086820446(9).(1) 

3. Results and analysis 

3.1. The lightcurve 

In Fig. [j] (upper panel) we show the mean light curve 
in the B-band measured with the MCCP at the 2.2m- 
telescope. In the lower panel of the diagram the phasing 
and phase coverage provided by our trailed spectra is dis¬ 
played. The most remarkable features of the light curve 
are: a) the deep eclipse of the white dwarf and the ac¬ 
cretion spot by the late-type companion star, b) a pro¬ 
nounced pre-eclipse dip centered on phase </> ec i= 0.88, c) 
a double-humped shape outside eclipse and dip, and d) 
strong QPO/flaring outside eclipse and dip. A detailed 
discussion of the light curves will be given elsewhere, in 
the context of our present analysis the interesting features 
are the eclipse and the dip. 


The eclipse is two-stepped with the first step being 
very steep (~3sec) followed by a more shallow second 
step with ingress time of about 120 sec. All four steps 
have different heights. We identify the obscured sources 
of light with the hot accretion plasma on the white dwarf 
emitting cyclotron radiation (steep ingress/egress) and the 
accretion stream (shallow ingress/egress). Both emitters 
radiate anisotropically which gives the explanation for 
the different stepsizes at ingress and egress. The eclipse 
of the white dwarf itself is not clearly resolved in our 
data. We measured the length of the eclipse as the dif¬ 
ference between the maxima in the first derivative of the 
lightcurve at ingress and egress phase to be 584.6 sec. 
This value represents the eclipse length of the accretion 
spot on the white dwarf surface. The spot is nearer to 
the secondary star than the center of mass of the white 
dwarf CM w( j. The eclipse of the CM w( j is shorter by a 
factor ocm/ ci^spot with a being the opening angle of the 
secondary star as seen from the CM and the spot, re¬ 
spectively. The minimum value of this correction factor is 
0.9792 (Atge^cM = 572.4sec), calculated for the extreme 
(and therefore unlikely) assumptions of inclination i = 90° 
and spot co-latitude (3 = 0° (spot in the orbital plane). 

The Gaussian shaped dip is centered on 4> ec \= 0.882 ± 
0.003 (deflection angle d = 42.5° ±1.1° with respect to the 
line joining both stars) and has a FWHM = 0.034 ± 0.001 
phase units. A ±2cr-range around dip center (95% of the 
absorbed flux) spans a range d = 32° — 53°. The emission 
lines are not obscured at dip phase (see Figs. ^ and [|) , only 
the cyclotron continuum is absorbed. The dip is therefore 
naturally explained as obscuration of the accretion spot 
by the transient accretion stream as in other polars (see 
Sect. 1). Since HU Aqr is a system with a high inclination, 
dip center and width mark exactly the location and extent 
of the coupling region in the orbital plane. 

3.2. The mean orbital spectrum 

In Fig. we show the mean-orbital, continuum- 
subtracted, radial-velocity-corrected spectrum of HU Aqr. 
For radial-velocity correction we used the sin-fit to the 
narrow emission line of He n A4686 (see below). The pro¬ 
files of the emission lines in Fig. || therefore show a spiky 
narrow emission line atop broad line wings. We can iden¬ 
tify the emission lines normally found in AM Herculis bi¬ 
naries, as e.g. the Balmer lines H/3 and Hy, the promi¬ 
nent lines of neutral (Hei AA4387, 4471, 4713, 4921, 5015) 
and ionized Helium (He n A4686), but find in addition 
several features not easily recognized in other systems. 
We note in particular the partially resolved Cm/Niil- 
structure at 4635-4650, Cn A4267, a second line of ionized 
He (HeiiA4541) with an additional unidentified fainter 
line in its red wing and MgnA4481 in the red wing of 
Hei A4471. 

The Ciil/Nm-complex is resolved into three lines, 
a single line of Nm (A4634.2) and two blends of Nm 
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Fig. 2. Mean-orbital high-resolution spectrum of HU Aqr on August 17, 1993, after subtraction of the continuum and after 
correction to zero velocity (narrow emission line). The same data are shown twice, in the bottom panel on an expanded intensity 
scale in order to facilitate identification of faint lines 


(AA4640.6,4641.9) and Cm (AA4647.4,4650.3,4651.5). 
The lines of C and N are of comparable strength. Ac¬ 
cording to Williams & Ferguson (1983) the Nm-lines are 
probably excited by continuum fluorescense. This does 
not apply to Chi A4650 or to Cli A4267 which do not 
couple directly to the ground states via electric dipole 
resonance lines. Whereas Cli A4267 may be excited by di- 
electronic recombination no final conclusion was reached 
by Williams & Ferguson about the probable excitation 
mechanism of the Cm-lines, apart from the general state¬ 
ment, that ‘some selective excitation mechanism is opera¬ 
tive for this line’. We would like to note simply here that 
this special excitation mechanism does not require other 
conditions regarding density, temperature and external ra¬ 
diation Held than are required for excitation of the other 
transitions of the high-ionization species found in the spec¬ 
trum of HU Aqr, since the kinematics and structure of all 
high-ionization lines seem to be the same. 

Three lines of the Heii-Pickering series (n = 4) lie 
in our spectral window at AA4859.32 (m = 8), 4541.59 
(9), and 4338.67 (10). The line at 454lA is clearly re¬ 


solved whereas the other two blend with H (3 and Hy. 
The wavelength-difference between the H-Balmer and the 
Hen-lines is 2.0lA and 1.80A, respectively, which is at 
the limit or below our spectral resolution at the speci¬ 
fied wavelengths. From our low-resolution spectra we esti¬ 
mate the flux ratio between the non- or weakly-blended 
lines at 4541.59A (m = 9) and 5411.52A (to = 7), 
F 5411 /F 4541 — 1.9, and interpolate to estimate the contri¬ 
bution of Hen A4859 to H /3, which is of the order of 10%, 
but variable as a function of phase. Hy itself is stronger 
than H (3 and the blending Hen-line (to = 10) is weaker 
than that at H/3, hence we expect a much weaker contam¬ 
ination of Hy. 

The Balmer lines show an inverted Balmer decrement 
with I(H(3) / 1(H y) = 0.86. The value for optically thin 
‘case B’ recombination is ~2.12 (Pottasch 1984), hence 
the Balmer lines are significantly optically thick. Although 
also HeiiA4686 is not optically thin (see below), the fact 
that this line is only weakly blended (by HeiA4713 and 
Cm A4650) and that our spectral resolution is highest just 
at this line makes it the most suitable one in our spectra 
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Fig. 3. Trailed, continuum-subtracted, high-resolution spectrum of HU Aqr on August 17, 1993, computed by phase-averaging 
the individual trailed spectra using 100 phase bins. The dynamic range of this grey-scale image with wavelength running along 
the abscissa and phase along the ordinate was set appropriately to emphasize weak lines. The behaviour of the strong lines is 
better visible on the zoomed representations of Fig. ^ 


for detailed investigations, but before doing so we take a 
general look at the trailed spectrograms. 

3.3. Trailed spectra 

In Figs. |] and [| we show the phased-averaged trailed 
spectra in overview (weak lines emphasized) and for the 
strong lines H/?, H 7 , Hei A4471, Hen A4686 in detail with 
a phase resolution of 100 phase bins. In Fig. || it can be 
seen that the phase intervals centered on <A> e d= 0.11, 0.46, 
and 0.78 which were covered only once during our obser¬ 
vations have higher noise in the continuum, the effect on 
the S/N in the strong lines is nevertheless negligible. All 
emission lines have the same structure and can be sub¬ 
divided into three emission line components (best visible 
in HeliA4686). Most striking is a narrow emission line 
(henceforth referred to as NEL), which is most prominent 
around (j> ec \= 0.5 and which has zero radial velocity at 
conjunction of the secondary star. This component is nat¬ 
urally explained as being of reprocessed origin from the 
illuminated hemisphere of the secondary star, seen here 
in unprecedented clearness. The second component has a 
more transient nature, it is a relatively narrow line, mov¬ 
ing in a zig-zag manner from maximum recessional veloc¬ 
ity (> 1000 kms -1 ) around eclipse phase to large negative 


velocities (< —1500kms -1 ) around ^ e d= 0.45. This com¬ 
ponent is brightest in the phase intervals <p e c i= 0.20 — 0.25 
and 0.70 — 0.75, around those phases where it crosses the 
NEL. We refer to this component as high-velocity com¬ 
ponent HVC. The third component is a broad underlying 
one visible throughout the whole orbital cycle (outside 
eclipse), with maximum bluesliift around </> e ci— 0.40 which 
we refer to as broad base component BBC. The HVC and 
the BBC must originate in the stream. 

Although the four main emission lines have the same 
general appearance differences exist in details. This con¬ 
cerns the width of the NEL, which is broader in the H- 
Balmer lines than in the He-lines, and the lightcurves of 
the lines (brightness variation as a function of phase). The 
width of the NEL of the different lines is addressed in 
the next section (3.4); the lightcurves are shown in Fig. || 
(integrated flux over the wavelength interval used for dis¬ 
play in Fig. [|). The intensity of the stream components 
dips to a minimum at phase </> ec i— 0.53, and is gener¬ 
ally stronger before than after this dip. This asymmetry 
is generally stronger for the low-ionization lines H 7 , H/3, 
HeiA4471 with flux ratios F(<j> < 0.53 )/F(<j) > 0.53) = 
2.09,2.14,2.02, respectively, than for the high-ionization 
He 11 A4686-line with the corresponding flux ratio of 1.54. 
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Fig. 5. Variation of the integrated line flux of the four main 
emission lines as a function of phase. Shown are from top to 
bottom the light curves for H 7 (histogram), H /3 (polygon), 
HeiiA4686 , and HciA4471. 

A likely explanation for the brightness variation of the 
HVC, which is pronounced only at certain phases, lies in 
the illumination geometry and significant optical thickness 
of the stream in the low-ionization lines. The Balrner lines 
and He 1 A4471 are bright when the stream is seen from 
its concave, that is its illuminated side (see Fig. [b| for a 
sketch of the geometry). The convex, non-illuminated side 
of the stream visible after </> ec i= 0.53 lies in its own shadow 
and is therefore much fainter than in the less optically 
thick Hell A4686-line. 

3-4- The width of the NEL in He- and H-lines 

The trailed spectrograms of the Balrner lines (in partic¬ 
ular their NELs) look somewhat blurred or diffuse when 
compared with the He-lines. Here we investigate, if this is 
caused by the weak blending lines of the Hell-Pickering se¬ 
ries mentioned above. The mean width (full width at half 
maximum FWHM) of the four strongest emission lines 
(H 7 , HeiA4471, HenA4686, H/3) in the phase interval 
= 0.50 — 0.56 (NELs intensive, stream components rela¬ 
tively faint) is (4.19±0.18)A, (2.74±0.15)A, (2.47±0.07)A, 
and (4.34±0.12)A respectively. The spectral resolution of 
arc lines at the wavelengths of our lines is 2.33A, 1.92A, 
1.67A, and 1.76A, which yields intrinsic widths of the NEL 
of 3.48A, 1.95A, 1.82A, and 3.97A. The He-lines have the 
same widths within the errors (1.95 vs. 1.82A), whereas 
the Balrner lines are significantly broader. 

We have synthesized blended line profiles of H j3 plus 
HeiiA4859 assuming the same intrinsic line width for both 
lines of 1.82A (that of the HeiiA4686 -line) and a maxi¬ 
mum relative intensity of 20% of the He-line. After folding 
with the instrumental resolution at H/3, the resulting un¬ 
resolved line had a FWHM of 3.8A (observed 4.34A). We 
conclude that the NELs of the Balrner lines are intrinsi¬ 


cally broader than those of the He-lines. Using the same 
scheme of synthesizing line profiles, we could estimate the 
width of the H/3-NEL to be 3.5A. Below we show that 
the width of the Hen-line is explained by velocity broad¬ 
ening introduced by the range of orbital velocities over 
the illuminated hemisphere of the secondary star. Hence, 
each line which originates in that region should have at 
least this width. If we assume that the observed excess 
width of H/3 is caused by thermal Doppler broadening, 
we would need temperatures in excess of 2 x 10 6 K in or¬ 
der to reflect the observations. This temperature, corre¬ 
sponding to coronal temperatures, is unbelievably high for 
a line originating in a (quasi-)chromosphere. The canon¬ 
ical temperature for the Balmer-line emission region in 
active M-stars is 10 4 K and the observed line widths at 
Ha are about ^ 2 A (Giampapa et al. 1982, Worden et 
al. 1981). This is also true for irradiated atmospheres of 
dMe stars although our knowledge of this very complex 
process is very poor (see e.g. Cram 1982). A more likely 
explanation seems to be that the observed line width re¬ 
flects orbital motion of a more extended Balrner line re¬ 
gion, perhaps some kind of wind/mass outflow component 
that co-rotates with the secondary star. An even simpler 
possibility is ‘saturation broadening’ due to large optical 
depths in the Balrner lines. 

3.5. Gaussian deconvolution of Hell \4686 

The conventional approach to investigate the emission line 
structure in polars is fitting of Gaussians. Due to the lack 
of spectral and temporal resolution often a distinction be¬ 
tween a broad and a narrow component is made and only 
two Gaussians are used but sometimes up to four (Rosen 
et al. 1987). We followed that approach and performed a 
triple-gaussian fit for the NEL, the HVC and the BBC in 
an iterative fashion. 

In a first step we fitted three Gaussians to each of 
the phase-binned spectra shown in Fig. [| while fixing the 
widths (FWHM) of the NEL and the HVC at 2.5Aand 
4.5A, respectively. This gave very good estimates for the 
radial velocities, and good estimates for the line inten¬ 
sities and the FWHM of the BBC. The latter displayed 
no distinct orbital variation, instead it showed an erratic 
scatter of 1.5-2 A around 11.5 A. Results for the radial ve¬ 
locities thus obtained are shown in Fig. The fits gave 
unreliable results for the line intensities and for the wave¬ 
lengths at those phases where two components are crossing 
each other. This happens at phases </> e d~0.23 and ^0.70 
between NEL and HVC and at </> e ci~0.35 between HVC 
and BBC. In order to provide input guess values for a sec¬ 
ond iteration, the phase-dependent radial velocities were 
approximated by smooth curves, a sinusoid for the NEL, 
spline-curves for the HVC and BBC. The sinusoid and the 
splines are shown as solid and dashed curves in Fig. |], the 
sine curve was calculated using only the data points drawn 
with filled symbols. In a second step we fitted all individ- 






Schwope, A.D. et al.: The eclipsing polar HU Aqr 


7 



Fig. 6. Radial velocities of the three subcomponents of 
He II A4686, determined by fitting Gaussians with fixed widths 
to the phase-binned spectra shown in Fig. ^ (NEL: open and 
filled circles, HVC: crosses, BBC: squares). The solid line is a 
sin-fit to the NEL-data shown with filled symbols, the dashed 
lines are smooth interpolations for the HVC and BBC used as 
guess-values for subsequent iterations 

ual 15 sec spectra in the trailed spectrograms, where these 
fits were restricted to a small range around the interpo¬ 
lated guess values and the NEL was fixed on the input 
guess values. Similarly, the phase-binned trailed spectro¬ 
gram was fitted a second time using these smoothed guess 
values from the first step but allowing the FWHM to vary 
freely within certain ranges. 

The results of these subsequent steps of the fitting pro¬ 
cedure, the intensity variations of the three emission line 
components and the variation of the FWHM of the NEL, 
are shown in Figs. |7| and In Fig. [Tojwe show only those 
data, for which the NEL could be clearly separated from 
the HVC. In Fig. |?] we also show the summed lightcurve 
of HVC and BBC representing the light from the stream. 
The lightcurve of the BBC can be roughly described as a 
two-stepped function, with constant brightness before and 
after superior conjunction of the secondary. It is brighter 
by ~50% during first half of the orbital cycle than dur¬ 
ing the second half. The hump centered on 4> ec \— 0.35 is 
an artefact of the fitting procedure, the excess light there 
is missing in the HVC-lightcurve, the same may apply to 
the less pronounced hump at phase 0.82. The lightcurve 
of the HVC has two pronounced humps separated by ~0.5 
in phase with the first hump centered on c^ ec i~0.22. The 
sum of both, the HVC and the BBC, shows two smooth 
humps. 

The lightcurve of the NEL finally is single-humped. 
This component is brightest around superior conjunction 
of the secondary, but has a remarkable asymmetry with 
respect to phase 0.5. It is brighter during first part of the 
orbital cycle, with maximum brightness at (f> ec i= 0.46. In 
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Fig. 7. Light curves of the three subcomponents of He II A4686, 
determined by Gaussian fits to individual trailed spectra as 
described in the text. The uppermost panel represents the 
summed light from the accretion stream, the lowest that from 
the secondary star. The solid line in the lowest panel is the 
result of a model calculation, described in Sect. 3.6 

some of our fitting experiments, depending on chosen in¬ 
put values and caging of parameters, the brightness of the 
NEL became even compatible with zero at phases later 
than 0.75 but it has ~20% of the peak flux at the sym¬ 
metric phase (j> e c i= 0.25. 

The width (FWHM) of the NEL is largest at and 
slightly after superior conjunction ( 4> ec \= 0.5), reaching a 
maximum value of about 2.45A. We were prevented from 
really measuring the width for phases later than ~0.63 due 
to the decreasing brightness of the NEL and interference 
with the other line components. 

The radial velocity variation of the NEL is in excel¬ 
lent agreement with a sin curve. The fit assuming a cir¬ 
cular orbit for the data points shown with filled sym- 
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bols in Fig. £j, v(cj>) = v sys + K 2 sin(27r(0 — 4>o)), yields 
v S ys = —7.8 ± 1.4kms _1 , K' 2 = 293.5 ± 1.3kms _1 , and 
(j )o = 0.016 ± 0.001. There is a small additional system¬ 
atic uncertainty not larger than ~3kms _1 for K 2 depen¬ 
dent on the data points included in the fit. If we cor¬ 
rect the systemic velocity to the solar system barycenter, 
Av = — 3.5 kms -1 , it becomes v sys = —11.3 ± 1.4kms . 
We refer to the radial velocity amplitude as K 2 (instead 
of I\ 2 ) because it represents the velocity of the center of 
light of the illuminated hemisphere instead of the center 
of mass. 

Glenn et al. (1994) have presented some 10 high- 
resolution spectra of Ha obtained when the system prob¬ 
ably was in a somewhat different state of activity. Proba¬ 
bly due to the lack of phase-resolution they discern be¬ 
tween a broad and a narrow component only and de¬ 
rive v sys = 83 ± 44 kms -1 , K = 244 ± 44 kms -1 , and 
cj )o = 0.04 ± 0.04, respectively, for their NEL. The results 
of both studies are compatible with each other only at the 
2 cr-level, hence we cannot completely rule out real changes 
of the emission line structure. Glenn et al. have described 
the broad component of Ha also with a sinusoid. Given 
our higher-quality data we cannot use such a simple ap¬ 
proach. 

Both of our stream components show maximum reces¬ 
sional velocity around eclipse, the full amplitude of the 
HVC is ~2000kms -1 , that of the BBC ~1400kms -1 . 
Neither the HVC nor the BBC display a radial velocity 
curve which can be described by a simple mathematical 
function. An analysis of these components requires a dif¬ 
ferent approach, which is presented below as tomogram 
analysis. If one, nevertheless, e.g. for comparison with the 
Glenn et al. analysis, fits a sinusoid to the radial veloc¬ 
ity data of the BBC, v(tft) = 7 + K sin(27r(<£ — 4>o)), 
one gets 7 = 3 ± lOkms -1 , K = 534 ± 15kms -1 , and 
</>o = —0.32 ± 0.01. Again, these values are compatible 
with the Glenn et al. analysis at the 1 2 a level. One ques¬ 
tion which is left open here in particular is, whether the 
two components HVC and BBC originating in the stream 
really reflect emission from two separate regions, as it ap¬ 
pears, or whether instead the emission all along the stream 
naturally produces skewed profiles simply due to optical 
depth, projection and curvature of the stream. The con¬ 
ventional approach of deconvolution using Gaussian fits, 
however, reveals very important information for the inter¬ 
pretation of the tomograms as far as e.g. intensity varia¬ 
tion is concerned. 

The asymmetrical shape of the lightcurve of the NEL, 
the non-zero velocity at eclipse and the non-vanishing 
phase offset <f>o are pronounced and require an interpre¬ 
tation. It appears likely that all three features are related 
to the same origin, an asymmetric illumination of the sec¬ 
ondary star. 


3.6. A reprocessing model for the NEL 

The skewness of the lightcurve of the NEL is significant. 
The sole explanation for such an asymmetry is to assume 
an asymmetric illumination of the secondary star. Sim¬ 
ilar asymmetries, although mirrored with respect to the 
axis of symmetry, have been observed by Beuermann & 
Thomas (1990) in the nova-like UX UMa-system IX Vel 
and by Hessman (1989) in the dwarf nova IP Peg. The 
likely mechanism for the asymmetry there is significant 
irradiation from the hot spot where the accretion stream 
impacts on the disk. Since this lies on the leading side of 
the secondary the trailing side is less irradiated and con¬ 
sequently fainter. Such an explanation cannot work in an 
AM Her binary, since the source of irradiation is bound to 
the surface of the white dwarf. 

A related type of asymmetry has been discovered in 
Nal absorption lines of AM Her itself by Southwell et 
al. (1995) and Davey & Smith (1996). These authors de¬ 
tect flux deficits on the trailing side of the secondary in 
the Nal absorption lines caused by stronger heating of the 
trailing hemisphere, and explain the necessary shielding 
of the leading side by absorption in the accretion column 
near the white dwarf, which points to the leading side in 
AM Her. 

We have synthesized line profiles originating from the 
illuminated hemisphere of the secondary star using the 
method as outlined in Horne & Schneider (1989) and 
Beuermann & Thomas (1990). In this model the Roche 
lobe surface is covered by a set of discrete surface ele¬ 
ments. The flux incident on each surface element is calcu¬ 
lated assuming a point source located at the position of the 
white dwarf. Lightcurves and line profiles are synthesized 
by summing contributions from all surface elements visi¬ 
ble at a given phase. At this stage we allowed for optically 
thin emission, being proportional to the incident flux, as 
well as for optically thick emission, where this flux was 
weighted with the projected surface of the element in con¬ 
cern. No further line broadening mechanism than velocity 
broadening over the illuminated Roche surface and the 
instrumental broadening were taken into account. We cal¬ 
culated a series of these profiles as a function of phase and 
shielding. We took the simple view of describing shielding 
by a weighting factor for each surface element of the lead¬ 
ing side which was constant for the whole hemisphere. We 
performed such calculations for three values of the mass 
ratio Q = M w d/M 2 = 3,4,5 assuming that the secondary 
follows the ZAMS mass-radius relationship given by Cail- 
lault & Patterson (1990) (in that case the mass of the sec¬ 
ondary is M 2 — 0.146Mq and only very weakly dependent 
on mass ratio) and for Q = 2.5, M 2 = 0.35 M 0 , as moti¬ 
vated by our tomogram analysis. The orbital inclination i 
was chosen according to the q — ^-relation for the observed 
length of the eclipse (see Fig. [If]) as i = 80°, 82°, 85°, 90°, 
respectively, for Q = 2.5,3,4, and 5. 
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Fig. 8. Results of model calculations for the NEL: Shown 
are the relation between radial velocity amplitude, systemic 
velocity and phase shift obtained from sine fits to synthetic 
NEL-profiles as a function of mass ratio and shielding of the 
leading hemisphere of the secondary star. The numbers des¬ 
ignate the mass ratio Q = The large open squares 

connected with solid lines are data computed in the optically 
thick mode, the small open squares connected with dashed lines 
are data computed in the optically thin mode (shown only for 
Q = 5) and the small circles are optically thin data for Q = 2.5 
with a massive secondary of M 2 = 0.35 Mg. Shielding increases 
by 10 % between the data points starting with no shielding at 
right. The large filled squares are observed data. 


The synthesized profiles were subsequently analyzed 
in the same manner as our observed spectra but first 
their spectral resolution was degraded to a FWHM = 
107kms _1 as observed. Then single Gaussians were fit¬ 
ted to determine the central wavelength. Afterwards sine 
curves were fitted to the modelled radial velocity data, 
taking into account only data from those phase intervals 
which were used also to determine the observed radial 
velocity curve. The results of this procedure are shown 
in Figs. H (parameters of radial velocity curves) and [iTi| 
(FWHM vs. phase) together with the observed data. The 
modelled radial velocity curves are not exact sinusoids, but 
the deviation is negligible for the phase intervals in con¬ 
cern. Two families of model curves are shown in Fig. |n| 
The upper four curves are calculated for Q = 2.5, M 2 = 
0.35 M@, the lower four for Q = 5,Af 2 = 0.146Mq (the 
other models with M 2 = 0.146 Mq but smaller Q are 
undistinguishable from that with Q = 5). The solid lines 
are optically thick models, the dashed and dotted lines op- 
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Fig. 9. Results of model calculations for the NEL: Shown are 
the relation between skewness (1(0 > O.5)//(0 < 0.5)) and 
phase of peak intensity of the NEL-lightcurve for a mass ratio 
Q = 5 and varaible shielding of the secondary star. As in Fig. |i], 
shielding grows from right to left in 10 %-steps and the larger 
symbols represent optically thick models. The filled square is 
the observed value with estimated uncertainties. 




Fig. 10. Full width at half maximum FWHM of the NEL 
transformed to velocity at the He 11 A4686 -line as a function 
of phase. The filled squares are observed data, the lines are 
model curves calculated assuming different mass ratios, mass 
of secondary star, optical thickness and shielding (see text for 
details). The dotted line at bottom indicates the instrumental 
resolution. 


tically thin. The upper curve of each pair is for no shield¬ 
ing, the lower curve for 60% shielding. 

In addition to the line profiles lightcurves for the repro¬ 
cessed lines as a function of phase were generated. They 
show a distinct asymmetry in the required manner, which 
was parametrized using (a) a skewness defined as ratio of 
fluxes F(cj) > 0.5 )/F{(j> < 0.5) and (b) the phase of maxi¬ 
mum intensity. Since these parameters are negligible func¬ 
tions of the mass ratio (they depend only on the adopted 
shielding), we show in Fig. [f]only one set of model data for 
Q = 5 and optically thin and thick modes, respectively, 
together with the parameters derived from the observed 
lightcurve of the NEL. 
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The models and data shown in Figs. || and || suggest 
consistently a mass ratio Q ~ 5, optically thick emission 
and an effective shielding of the leading side of the sec¬ 
ondary star by 60%. Lightcurves for the optically thick 
and thin models were shown also by Horne & Schneider 
(1989), where it can be seen that the thin models predict 
a plateau around </> ec i= 0.5. This happens when always 
the same numbers of surface elements are visible for a dis¬ 
tinct phase interval. Such a plateau is not observed in HU 
Aqr and we conclude from all the detailed comparisons 
between our models and the observed data, that the NEL 
is non-negligibly optically thick. 

The fit is not completely satisfactory for the variation 
of the FWHM. While the low-Q, high-M 2 model predicts 
too high values for the FWHM for most orbital phases 
(upper four curves) the model with Q = 5 predicts some¬ 
what too small values, at least if one takes shielding into 
account and uses the optically thick models. In addition, 
shielding moves the phase of maximum FWHM from (f> ec i= 
0.5 towards smaller values whereas the observed FWHM 
peaks only at </> e d— 0.5 — 0.55. We note, however, that 
the FWHM is the least well determined quantity among 
all those used here for comparison between observation 
and models and we do not regard the deviation between 
observation and simulation as a severe problem. 

3. 7. Tomogram analysis of the four main emission lines 

In this section we investigate the trailed spectrograms 
shown in Fig. [| using a tomographic inversion technique 
known as filtered backprojection and described e.g. by 
Robinson et al. (1993). Making use of such an inver¬ 
sion technique implies that any emission not observed 
at rest wavelength of a line is regarded as shifted due 
to the Doppler effect and the backprojection from the 
observers (v, ((>)-space to the original (v x , t> y )-plane thus 
makes use of the whole line profile (instead making use of 
only e.g. width, intensity or center). The interpretation of 
such maps is mostly done by comparison with computer¬ 
generated schematic drawings of the possible locations of 
line emission in the binary using as input parameters only 
the velocites of the two stars of the binary. 

The computation of a predicative Doppler-tomogram 
(or Doppler-map or shortly map) requires a large number 
of projections (high phase resolution) and a high spectral 
resolution. Our dataset, in particular the trailed spectro¬ 
gram of Hen A4686, is unique in this respect. The tomo¬ 
grams shown as grey-scale images in Fig. [TT| and as con- 
tourplots in Fig. [1^ have a velocity sampling of 5kms _1 
for HenA4686 and 20kms _1 for the other lines. Two 
features in the maps are striking, a bright spot at veloci¬ 
ties v x ~ 0 — 30kms _1 , v y ~ 290kms _1 and a cometary 
tail linked to it and extending to v x ^ — lOOOkms -1 at 
v y — Okms . A third feature appears as faint smeared 
blotch in the lower left quadrant of the maps. Since this 
blotch has a low contrast in Fig. |ll] only we show in 


Fig. |l2] the Doppler map of those emission components 
of He IIA4686 attributed to the accretion stream. It was 
computed by subtracting the Doppler image of the NEL 
from that of the complete profile and shows two main fea¬ 
tures, the cometary tail and the blotch below it. 

Backprojection transforms sine curves to points, the 
bright spot thus corresponds to the NEL and repre¬ 
sents the irradiated secondary star. For the He 11 A4686- 
line we compare the position of the spot in the Doppler 
map with the results of our Gaussian fits. The spot 
appears at v x = 30kms _1 and v y = 283kms _1 . The 
offset of the spot with respect to the symmetry axis 
v x = 0.0 reflects the obscuration of the leading (left) 
side. The ratio of the velocity components corresponds 
to the measured phase offset </>o )S m of the radial veloc¬ 
ity curve (f> o,dop = (27r) -1 axctan{v x / v y ) = 0.017, in good 
agreement with </>o,sin = —0.016. The velocity amplitude 

^ 2 ,do P = \J v x + v y = 285, however, is lower than found 
from the sine fit to the radial velocity curve K' 2 sin = 
293.5 km s _1 . This might be caused by deviations of the 
true line profiles from Gaussian curves. We see such de¬ 
viations pronounced in our modelled line profiles in the 
optically thin approximation and they are present too but 
less pronounced in the optically thick models (for an il¬ 
lustration see Horne & Schneider, 1989, their Fig. 12). We 
have checked the hypothesis by comparing radial veloci¬ 
ties measured (a) by sine fitting radial velocity curves and 
(b) by 2D-Gaussian fitting in tomograms of our modelled 
profiles. We indeed find differences of the required type 
of the order 3-8kms _1 depending on mass ratio and opti¬ 
cally thin/thick approximation. An additional small bias 
towards lower velocities of the observed peak in the tomo¬ 
gram is likely due to the underlying cometary tail. So we 
regard both methods to be in agreement. The maps cal¬ 
culated for our synthesized NEL-profiles show that even 
with our low spectral resolution (with respect to the width 
of the NEL which is not well sampled) we should have seen 
a triangle-like structure in the Doppler map sketching the 
heated ‘nose’ of the secondary star if the emission hap¬ 
pens in the optically thin approximation. In the optically 
thick approximation the modelled tomogram yields only 
a slightly elongated structure, very similar to what is ob¬ 
served. 

The Roche lobes shown as overlays in Fig. [□] have all 
the same size thus facilitating a comparison of the four 
emission lines. The Roche lobes were calculated for a mass 
ratio of Q = M w< j/M 2 = 4. According to the larger widths 
of the NELs of the Balmer lines discussed in Sect. 3.4 the 
corresponding spots are more smeared than those of the 
Helium lines. In addition, emission is more concentrated 
towards the equator of the Roche lobe in the Balmer lines. 
All spots except that in the map of H/3 are bright on 
the trailing side. The H/3-spot has maximum intensity 
at v x — O.Okms -1 and a longish extension inclined by 
135° towards negative (v x , n y )-velocities. The radial veloc- 
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Fig. 12. Doppler map of the stream components of He II A4686 
computed by sutracting the Doppler image of the NEL from the 
Doppler image of the whole line profile. In order to emphasize 
the blotch at (v x ,v y ) ~ (—400, —300) kms -1 , we adopted a 
linear intensity scale between 5% and 60% of the peak intensity. 
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ity amplitudes of the NELs of HeiA4471, H 7 and H/3 as 
measured in the tomograms are v = 284,313, 312kms _1 , 
respectively. The value for He 1 A4471 agrees with that of 
HellA4686, both values for the Balmer lines agree with 
each other but are significantly higher than those for the 
Helium lines. 

The cometary tail connected to the secondary star in 
the tomogram is the horizontal accretion stream. It is the 
Doppler image of the HVC as can be seen directly by 
mapping only this component singled out using the Gaus¬ 
sian fits. Again this feature appears most brilliant in the 
map of Hell A4686. Its location and curvature in the map 
are sensitive to the velocities of the stars in the binary 
and, hence, to the masses and below we will give mass es¬ 
timates relying on the assumption that the ridge of light 
exactly reflects the location of the trajectory of the stream 
(approximated in the usual manner as single particle tra¬ 
jectory, Lubow & Shu 1975, freely falling in the gravitional 
potential). 

The third component in the Doppler-map, the diffuse 
blotch in the lower left quadrant, can be identified with 
the BBC in the trailed spectrogram. This again became 
evident by mapping the BBC alone which was isolated us¬ 
ing the Gaussian fits (the corresponding center of emission 
of the BBC is marked by a rhombus in Figs, [fl] and |l3| ). 
The blotch represents emission from the vertical stream. 

We have modelled the location of all components in the 
tomogram in a manner described e.g. by Marsh & Horne 
(1988) but included in addition a magnetic drag influenc¬ 
ing the trajectory of the stream. The results are shown for 


two extreme cases of the mass ratio in Fig. [k] (Q = 2.5 
and 5, respectively) and for an intermediate value of Q = 4 
in Fig. |Il|. The corresponding picture in positional coordi¬ 
nates is shown in Fig. |l4|. In Fig. [l^ we show the geometry 
of the stream as seen by a hypothetic observer at the white 
dwarf. 

The shape of the secondary star projected onto the 
orbital plane is the same in position and Doppler coordi¬ 
nates, though rotated by 90° in the Doppler coordinates 
because the orbital velocity vectors are perpendicular to 
the position vector. Our stream starts at the inner La- 
grangian point L\ with no initial width. The stream then 
follows a single particle trajectory (Lubow & Shu 1975) 
in the Roche-potential down to the interaction region be¬ 
tween magnetosphere and accretion stream. The location 
and size of this stagnation region is determined from the 
phase and the phase width of the absorption dip (Fig. ^]). 
In Figs. H to |l5| we show three trajectories with turn¬ 
off points adjusted to reflect the center and width of the 
stream (distance to center is ±2cr of Gaussian fitted dip). 
There the stream leaves the orbital plane as it couples 
onto field lines. We assumed no dissipation to occur in 
the interaction region. The abrupt re-orientation of the 
stream happens in a small spatial region, but produces a 
large jump in velocity from ~(-1100,0) to ~(-400,-300) in 
(v x ,v y )- space. The stream then follows the field lines of 
an assumed dipolar field. The field axis was fixed at colat¬ 
itude (3 = 25° measured with respect to the rotation axis 
and azimuth 'F = 40° with respect to the line joining both 
stars. The former value derives from the observed motion 
of cyclotron lines in low-resolution spectra (Schwope et 
al. 1996, in preparation), the latter one from the phase of 
the optical/X-ray bright phase center. 

We have chosen a model with Q = 4 (K\ = 87 kms -1 , 
K 2 = 350 kins -1 ) for the overlay in Fig. |ll] which gives 
a fairly good representation of the four lines and their 
subcomponents regarded there. We note, however, that 
the HVC-emission of the HenA4686-, HeiA4471-, and 
H/3-lines is concentrated below the line defined by the 
model, whereas that of the Hy-line is concentrated above 
the model line. We find excellent agreement between the 
HVC-ridge of HeilA4686 and the trajectory for a mass 
ratio of Q = 2.5 (A'i = 147kms -1 , K 2 = 367kms -1 ) 
and good agreement for H 7 using a mass ratio of Q = 5 
(Ad = 77 kms -1 , K 2 = 384kms -1 ) (Fig.|I|). 

We see in the tomogram that significant broadening 
takes place all along the horizontal stream, which is not 
reflected by our model. This is due to the simplicity of 
the model which does not take into account any poten¬ 
tial broadening mechanism as e.g. the finite width of the 
stream and its optical thickness. 

It is somewhat surprising that we cannot detect 
emission from the stream at really high velocities (> 
700 kms -1 ) close to the white dwarf (ve — 4000 kms -1 ). 
This may reflect the effects of a high degree of ionization 
and small radiating surface of the stream at high veloci- 
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Fig. 14. Scheme of the accretion geometry in HU Aqr pro¬ 
jected onto the orbital plane for an assumed mass ratio Q = 5. 
Shown are the loci of the center of mass (cross), the white 
dwarf (circle) and the secondary star (cross at right) as well as 
the Roche-lobes of both stars. Matter leaves the secondary star 
at the Li—point and follows a ballistic trajectory. Three held 
lines connecting the stagnation region and the white dwarf are 
shown, the region in between is density dotted. One additional 
held line (dashed line) connects the inner Lagrangian point 
with the white dwarf. The region between this line and the 
dense vertical stream, where the accretion curtain is located, 
is dotted lightly. 


Fig. 15. Scheme of the accretion geometry in HU Aqr as seen 
by a hypothetical observer located at the center of the white 
dwarf. Matter leaves the secondary star at the Li—point and 
follows a ballistic trajectory, hence it is dehected under the 
influence of Coriolis forces to the left. The stream leaves the 
orbital plane and couples onto held lines at an azimuth of 43°, 
where the absorption dip occurs. Three held lines connecting 
the stagnation region and the white dwarf are shown indicat¬ 
ing the width of the absorption dip (dark dotted region). One 
additional trajectory (dashed line) is shown leaving the plane 
at the Li—point. The region between this line and the dense 
vertical stream is less dense dotted (accretion curtain). 


ties. Possible effects which in addition may contribute to 
the low velocities in the BBC (and its Doppler image, the 
blotch) are (a) a slightly incorrect continuum subtraction 
thus removing faint emission at high velocities, and (b) 
emission not only from the accretion stream but also from 
the accretion curtain, hence from matter leaving the plane 
all along the horizontal stream. This matter is thought to 
be responsible for shadowing of the leading side of the 
secondary star and its re-emission has not been regarded 
so far. The turn-off points of such trajectories lie closer 
to the L\ and the jumps from the HVC- to the BBC- 
region in the map are smaller thus closing the expected 
gap. One sees also that the observed large width of the 
BBC is introduced naturally by the models due to the sim¬ 
ple fact that coupling happens along the trajectory over a 
significant distance and the jumps from the HVC- to the 
BBC-region adjust correspondingly. In principle, so-called 
‘gamma-smearing’ may contribute to the observed width 
of the BBC in the Doppler map (Diaz & Steiner 1994). 
Gas with velocity perpendicular to the orbital plane v z 
produces a shift 7 = v~ cos i. The effect of this is to blur 
the corresponding feature in the Doppler map over a ring 
of radius 7 centered at its usual position. If we assume 
v z = 1000kms _1 and i in the range 81°-85°, we get a 
smearing amplitude of 156 — 87kms _1 , which is of the or¬ 
der of or below the resolution of our map. Hence, our map 
of HU Aqr is essentially blind to the vertical motions. 


In Figs. 14 and [U4 the distribution of matter in the 


system is shown. The dark shaded region in both figures 
are the horizontal and vertical streams, while the dotted 
region gives a schematic of the accretion curtain neces¬ 
sary for shielding the leading hemisphere. This curtain is 
constrained by a trajectory which couples to the magnetic 
held at the Li-point. This gives not the required 60% de¬ 
rived above obscuration indicating that the view taken 
here is still too simplistic. Additional obscuration will be 
provided by nature by the finite width of the stream at 
the L\ which is of the order of 10% of the secondary star’s 
radius. 


We can use our stream models to obtain estimates for 
physical parameters of the stagnation region. The extent 
of this region in the orbital plane (length along the tra¬ 
jectory) is l ~ 6 x 10 ° cm and along which the undis¬ 
turbed dipole held strength varies between 3 and 9 kG. 
The polar held strength was set to B = 38 MG which 
yields a held strength in the accretion spot of about 36- 
37 MG as required by the identification of cyclotron har¬ 
monics (Schwope et al. 1993, Glenn et al. 1994). Our code 
assumes a constant cross section along the trajectory, 
pv = rh = constant. Coupling at different loci along the 
trajectory is possible by a variation of the density p , thus 
varying the ram pressure p ram oc pv 2 of the stream. Un¬ 
der these assumptions we found that the ratio of the ram 
and magnetic pressure varies by a factor of 6-7 along the 
trajectory in the stagnation region. The size of the stagna- 
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tion region is transferred down to the white dwarf surface 
via dipole field lines where it occupies a region with full 
opening angle of about 5°-6°corresponding to a linear di¬ 
mension of l ~ 7 x 10 7 cm for a white dwarf with radius 
-Rwd = 7 x 10 s cm. 

The tomogram together with the schemes of the ac¬ 
cretion geometry also yields hints for an understanding of 
the lightcurves of the HVC and the BBC (Fig. [?]). Both 
components are brighter during first half of the orbital 
cycle, when we view the magnetic curtain’s concave il¬ 
luminated side. The HVC is bright and easily to recog¬ 
nize in the trailed spectrogram around quadrature phases. 
The observed spectra at these phases correspond to pro¬ 
jections of the map parallel to the v x — axis, hence along 
the cometary tail, producing a sharp feature in a spectro¬ 
gram. The observed spectra at conjunction phases, on the 
other hand, correspond to projections along the v y — axis. 
At these phases emission from the horizontal stream is 
smeared over an interval of ~1000kms -1 which produces 
only very broad spectral features. Hence, at theses phases 
the HVC apparently vanishes and we obtain the minima 
in the lightcurve of the HVC. 


4. Discussion and conclusions 

We have found two distinct features in the trailed spectro¬ 
grams of the four main emission lines in the blue spectral 
regime of HU Aqr, in particular the line of ionized Helium 
HellA4686, which are potentially useful for a determina¬ 
tion of the stellar masses. These are the amplitude of the 
radial velocity curve of the narrow emission line NEL and 
the location and curvature of the horizontal stream in the 
tomogram, the Doppler image of the HVC. An additional 
constraint on the mass ratio is set by the observed eclipse 
length. All these constraints are shown in a (Q,i)— plane 
in Fig. |l(] with Q being the mass ratio and i the orbital 
inclination. The eclipse length is dependent on the size of 
the secondary star’s Roche lobe and the inclination, hence 
relation (1) in the diagram is purely geometrical (it is of 
course implicitly assumed that the secondary fills exactly 
its Roche lobe) (Chanan et al. 1976). The solid line is valid 
for the measured eclipse length, the nearby dashed line is 
valid for the maximum correction which accounts for the 
offset of the accretion spot from the white dwarfs center, 
see above. Relation (2) in the diagram for the NEL makes 
use of a factor which relates the center of light of the 
illuminated hemisphere to the center of mass. This fac¬ 
tor is dependent on Q and shown in Horne & Schneider 
(1989) for the optically thin and in Schwope et al. (1993) 
for the optically thick models. The differences between 
both modes are negligible as far as this factor is concerned. 
The computation of the curves (2) in the (Q,i )—plane 
makes use further of the Roche geometry (approximation 
of the spherical equivalent Roche-radius as a function of 
Q given by Eggleton (1983)) and Kepler’s third law. We 
have assumed that the secondary is a main sequence star 



Fig. 16. Mass determination of HU Aqr using the observed 
length of the eclipse (1), the radial velocity amplitude of the 
illuminated hemisphere of the secondary star with different 
mass-radius relations (2), and the curvature of the ballistic 
stream in the Doppler map (3). 


and can be described by a ZAMS mass-radius relation for 
late-type stars. With M 2 = 0.146 Mq, the empirical re¬ 
lation by Caillault & Patterson (1990) predicts the least 
compact secondary for HU Aqr, the other extreme with 
M 2 = 0.206 M 0 is due to the relation by VandenBerg 
et al. (1983) (solid and dashed lines in Fig. [l(]). Relation 
(3) finally is the result of a fit-per-eye to the location and 
curvature of the horizontal stream in the tomogram which 
relies on the assumption that the center of light along the 
ridge of HeiiA4686 marks the single-particle trajectory. 
This results in the smallest estimates of Q with Q ~ 2.5. 


It is clear from Fig. 16 (as well as from several pre¬ 
vious Figures) that no unique solution exists presently in 
the framework of the models presented here. Relations (1) 
and (2) suggest high values of the mass ratio and the in¬ 
clination with Q ~ 4.5 — 5.5, i > 85° for the Helium lines 
and Q = 5 — 5.5,* > 87° for the Balmer lines (the latter 
are not shown in Fig. [ui]). The curvature of the stream 
of Hen A4686, HeiA4471 and H (3 on the other hand sug¬ 
gests, that the mass ratio might be as low as Q = 2.5 and 
that the inclination correspondingly might be as low as 
i ~ 80°. Both estimates of the mass ratio seem to exclude 
each other presently and it appears likely that neither 
the NEL (in combination with AVcorrection and ZAMS- 
approximation) nor the HVC (fitted by a single-particle 
trajectory) can be used straightly to give an accurate mass 
estimate. 
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The disagreement between both approaches becomes 
smaller if one assumes that NEL emission is more concen¬ 
trated towards the center of mass of the secondary than 
we have assumed. An absorbing cloud of matter in the 
Li-region which would suppress emission from the ‘nose’ 
of the secondary could e.g. provide the necessary bias of 
the effective radial velocity amplitude towards higher val¬ 
ues (and could also lower the width of the line profiles, 
Fig. |Ty). The expected re-emission might be responsible 
for the longish extension of the bright spot near the L\- 
point in the tomogram of H (3. That such bias exists ap¬ 
pears likely from the observed higher NEL-velocity of the 
Hydrogen Balmer lines with respect to the Helium lines. 
The ‘nose’ is more effectively shielded in the more opti¬ 
cally thick Hydrogen lines. 

The best way to solve the conflict and to calibrate the 
M/R-relation of the secondary star would be to measure 
the radial velocity of the secondary star independently 
using e.g. photospheric absorption lines. This will help to 
sort out which feature (NEL, HVC) of which line may 
serve as the estimator of the mass ratio. It requires HU Aqr 
to be observed in a low state of accretion. Such a measure¬ 
ment is highly demanded in order to calibrate either of the 
methods used here for mass determination. Another hint 
to the likely mass of the stars and the mass ratio could 
be derived from ingress/egress length into/from eclipse of 
the white dwarf itself, rather than ingress/egress of the 
accretion spot on it. 

Our model with Q = 2.5 implies that the secondary 
star in HU Aqr must be as massive as 0.35Mg. It would 
be much more compact than a ZAMS-star and in con¬ 
sequence much more luminous. This would imply that 
the distance estimates given by Glenn et al. (1994) and 
Schwope et al. (1993) are lower limits only. 

We have presented high-resolution spectral observa¬ 
tions of the eclipsing AM Herculis star HU Aqr obtained 
when the system was in a high state of accretion. We have 
identified three emission line components in our trailed 
spectra of He n A4686, He IA4471, and the H-Balmer lines. 
These could be uniquely identified as originating on the 
illuminated hemisphere of the secondary star (the narrow 
emission line NEL), the stream in the orbital plane (the 
high-velocity component HVC) following a ballistic tra¬ 
jectory, and the stream coupled onto field lines out of the 
orbital plane (the broad-base component BBC). We have 
uniquely identified the horizontal stream in this AM Her¬ 
binary by means of Doppler tomography and could trace 
it from the Li-point down to the stagnation region. On 
the assumption that the NEL of He n A4686 is formed on 
the surface of the Roche-lobe of the companion star we 
have determined its mass. This estimate is in disagree¬ 
ment with an estimate derived from the tomogram. We 
propose that Doppler tomography of the ballistic stream 
can be used as a new tool for the mass determination of 


accreting binaries provided we can calibrate this method 
using an independent mass determination. 
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Fig. 4. Grey-scale representation of trailed, continuum-subtracted, high-resolution spectra of the strong emission lines of HU Aqr 
on August 17, 1993 (zoomed representation of Fig. ]|). Phase runs along the ordinate from bottom to top, wavelength has been 
transformed to velocity using the rest wavelengths of the specified lines. The dynamic range for all representations is for all 
lines approximately between 0 and 80% of maximum intensity. 
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Fig. 11. Doppler maps of the four main emission lines of HU Aqr computed by filtered backprojection of the trailed spectra 
shown in Fig. The velocity resolution for the H 7 -, Hel-, Hell-, and H/3-lines is 161, 129, 107, and 109kms -1 , respectively, 
at the origin and 173, 143, 124, 126kms” 1 at v = 1000kms -1 . The overlay shows the shapes of the secondary star and the 
accretion stream in Doppler coordinates for an assumed mass ratio of Q = M w( j/A 72 = 4. 
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Fig. 13. Contour plots of the Doppler maps for HellA4686 (left) and H 7 (right) shown in Fig. Ill with 11 equidistant contours 
between 0% and 100% of the peak intensity. The rhombus at ( v x ,v y ) ~ (—400, —300) km s~~indicates the center of BBC 
emission, the small circle and the cross give the loci of the centers of mass of the white dwarf and the secondary star for the 
assumed mass ratios of Q = 2.5 (left diagram) and 5 (right diagram). Shown are the shape of the secondary star, the ballistic 
trajectory (horizontal stream, upper left quadrant) and the part of the stream which has left the orbital plane (vertical stream, 
lower left quadrant). 










